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CN ; I. INTRODUCTION 

, ^ , Compton scattering, 7*+p — >■ 7+p, with initial photons real or virtual, has been intensively investigated theoretically 
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^> ' and experimentally |1|-|17|. In the case of deeply- virtual Compton scattering (DVCS), where the initial photon is 

' highly- virtual, the QCD factorization has been proven [5, 8, 9] and the amplitude can be expressed in terms of the 

5~, ■ generalized parton distributions (GPD) [1-^, 8-16] convoluted with some hard coefficient function. However, it is not 

' possible to make a deconvolution and unambiguously extract the GPDs of the target from the experimental DVCS 

1—1, data. The standard procedure in this case is to construct a plausible model with some free parameters, fix the free 

i-G ' parameters from the experimental data and after that make conclusions about the GPDs of the target. 

^p-^i The situation is different in case of double deeply virtual Compton scattering (DDVCS), ^* +p — >■ 7*-t-p — >■ ll+p \\M - 

O ■ |20 |. As it has been shown in [T^, the DDVCS amplitude allows such a deconvolution. Unfortunately, the corresponding 

D \ cross-section of the process is suppressed by aem/37r compared to ordinary DVCS and falls into a picobarn- level range, 

i-G ' so it cannot be accessed by the existing accelerating facilities. However, at future accelerating facilities [2l|, [22| with 

higher luminosities, this cross-section can be measured. A special case, when the initial photon is real, has been 

fT^ . studied in [2J, [2J| , and in the color dipole framework in [2^ . This process may be viewed as time- inverted DVCS with 

J> ■ negative virtuality —Mj, so the color dipole model is applicable for high lepton masses, Mj ^ M'^. A high luminosity 

^~~^ , flux of real photons may be created by ultraperipheral scattering of protons on nuclei; such collisions at the LHC 

CN ■ will be an ideal tool for the study of this process. However, it is not possible to access the whole DDVCS amplitude 

in ultraperipheral hadronic collisions, since the virtuality of the initial photon, controlled by the formfactors of the 

colliding hadrons, cannot be very high. Thus, the only way to study the DDVCS is the electroproduction of lepton 

'O ' pairs, e+p— > e + p + II. 

Recently DVCS on proton and nuclear targets has been studied within the color dipole approach in |26l - [3(j |. In 
this paper we extend that study to the DDVCS case. We focus on the high-energy kinematics range, since the 
future accelerating facilities are expected to operate at these energies [2l|, [22], where the effects of coherence are 
important. The general framework for evaluating the shadowing corrections is the Gribov- Glauber approach |3l|,|32|. 
In particular, we take into account gluon shadowing corrections, which correspond to the triple-Pomeron diffraction 
in the Gribov inelastic shadowing corrections. These corrections onset at xb ^ 10~^, and give a sizeable contribution 
C^ ^ to nuclear shadowing at xb ~ 10"*^. 

The paper is organized as follows. In Sections HIl we review the general formalism of the color dipole approach. In 
Section IlIII we discuss the frequently used frozen approximation, which is valid for asymptotically large energies. In 
Section IIVI we discuss the method which will be used for calculations of nuclear shadowing effects and demonstrate 
that for asymptotically large energies it reproduces the results from Section IIIIl In Section |V] we discuss the gluon 
shadowing and its effect on the DDVCS observables. In Section IVTl the wavefunction of a real photon is evaluated in 
the instanton vacuum model. In Section [VIII we present the results of numerical evaluation and draw conclusions. 
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II. COLOR DIPOLE MODEL 

The electroproduction cross-section for unpolarized lepton pairs in the DDVCS process has a form 



dtdMl dt SttM? \ M? 



where da^f /dt is the cross-section of the process 7*p — > 7 *p, Mj. is the invariant mass of the produced lepton pair, 
and t is the square of the momentum transferred to a target. In what fohows, we wiU concentrate on evaluation of 
the cross-section dajf /dt. For the unpolarized Compton scattering discussed in this paper we obtain, 
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where the amplitude Ajll is evaluated in the color dipole model. In this model the dominant contribution to the 
Compton amplitude comes from gluonic exchanges. Then the general expression for the Compton amplitude on a 
nucleon has the form, 

1 
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A^P (s. A, Q^,Aq) = eWeW ^ d/i^d^^d^v^d\^^f [fi^, ^2) A'' (/3i, rl; Z^^, r!,; A) vI/W (^p^^^^) , (3) 



where e)l is the photon polarization vector; /3i_2 are the light-cone fractional momenta of the quark and antiquark, ri 2 
are the transverse distances in the final and initial dipoles respectively; A is the momentum transfer in the Compton 
scattering, A'^{...) is the scattering amplitude of the dipole on the target (proton or nucleus), and ^^.^, ,. {/3, r) are the 
light-cone distribution functions of the initial and final photons in the polarization state i . When the virtuality of 
the corresponding photon is large, we may use the well-known pQED expressions [33|, |3Jj, 

^f *m - ((«' + (1 - «)') eie2i^i(eir)X2(e2r) + mlKoieir)Ko{e2r)) , 



where e^ — ^ Q^a{l — a) + rn^, and make an analytical continuation of the space-like wave functions to the time-like 
photons by a simple extrapolation, Q^ -^ ^^fi 123; ^^'^ ^^^ 

Ko{ix) = --{Ux)-iY^{x)), (4) 

K^{ix) = ~^{Ji{x)-iY^{x)), (5) 



where x = r ^ Mj a{l — a) — rrvi for the kinematics Mj. 3> m'i/{a{\ — a)). However when the virtualities are small, 

we have to resort to some model (see Section I VII for more details) . 

At high energies in the small angle approximation, A/y^ ^ 1, the quark separation and fractional momenta /3 are 
preserved, so 
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is the gluon GPD of the target, P' = P + A, P = (P + P')/2, G^,y{x) is the ghion loop operator, £00 {x, y) is the 
Wilson factor required by gauge covariance. For this GPD we use a gaussian parameterization (Sa - ISq . 
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where the phenomenological functions (To(a;), Pg(a:), i?(a;) are fitted to DIS, real photoproduction and i:p scattering 
data. We discuss them in more detail in Section IVlIl The parameterization ^ does not depend on the longitudinal 
momentum transfer and decreases exponentially as a function of A^. Since the parameterization (|9]) is an effective one 
and is valid only in the small- a; region, we do not assume that it satisfies general requirements, such as positivity [39| 
and polynomiality [3] constraints. 

The prefactor {k + y) (fc — ■§■) in © guarantees convergence of the integrals in the parameterization ^. In the 
forward limit, A — ^ 0, the amplitude © reduces to the saturated parameterization of the dipole amplitude proposed 
by Golec-Biernat and Wlisthoff (GBW) |40|, 



ad(/3,r) = 2Jd%'^mfg{r,b',P)^^j'^a,{k^)F(xXo) 

I ampli 
by pQCD |4l| 
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Generally, the amplitude /i^(...) involves nonperturbative physics, but its asymptotic behavior for small r is controlled 



/^;(r,A,/3),^o(xr^ 

up to slowly varying corrections ~ ln(r). 

The calculation of the differential cross section also involves the real part of scattering amplitude, whose relation 

to the imaginary part is quite straightforward. According to [42|, if the limit lim ( ™-' ) exists and is finite, then 

the real and imaginary parts of the forward amplitude are related as 



7^e/(A = 0) = s"tan 
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In the model under consideration the imaginary part of the forward dipole amplitude indeed has a power dependence 
on energy, Imf{A = 0; s) ~ s", so ([TT|) simplifies to 
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= tan(f(a- 1)) = e. 



(12) 



This fixes the phase of the forward Compton amplitude, which we retain for nonzero momentum transfers, assuming 
similar dependences for the real and imaginary parts. Finally we arrive at. 



yW^ 



,U), 



A'^^l^ = {e + t)e^Hq')ey\q) / d'r J dp-^^ {fi,r)^Ynifi.r)^^ f^^{r, \ P,s), 





(13) 



III. NUCLEAR SHADOWING IN THE FROZEN LIMIT 



Nuclear shadowing signals the closeness of the unitarity limit. Hard reactions possess this feature only if they have 
a contribution from soft interactions. In DIS and DVCS the soft contribution arises from the so called aligned jet 
configurations [43|, corresponding to qq fluctuations very asymmetric in sharing the photon momentum, /? ^ 1. Such 
virtual photon fluctuations, having large transverse separation, are the source of shadowing |49| . 

Calculation of nuclear shadowing simplifies considerably in the case of long coherence length [4J], i.e. long lifetime of 
the photon fluctuations, when it considerably exceeds the nuclear size. In this case Lorentz time dilation "freezes" the 
transverse size of the fluctuation during propagation though the nucleus. Then the Compton amplitude of coherent 
scattering, which leaves the nucleus intact, has the same form as Eq. (|13l) with a replacement of the nucleon Compton 
amplitude by the nuclear one. 
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where b is impact parameter of the photon- nucleus colhsion, Ta(&) = /_ dz PA{b, z) is the nuclear thickness function, 
given by the integral of nuclear density along the direction of the collisions. In this expression we neglect the real 
part of the amplitude which is particularly small for a coherent nuclear interaction. 

For incoherent Compton scattering, which results in nuclear fragmentation without particle production (quasielastic 
scattering), the cross section has the form pa |. 
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Here i?e; is the i-slope of elastic dipole-nucleon amplitude. In this equation we treated the term quadratic in the 
dipole cross section as a small number and expanded the exponential in curly brackets. 

IV. ONSET OF NUCLEAR SHADOWING 

The regime of frozen dipole size discussed in the previous section is valid only at very small xb, or at high energies. 
However, at medium small xb a dipole can "breath", i.e. vary its size, during propagation through the nucleus, and 
one should rely on a more sophisticated approach. 

In this paper we employ the description of the onset of shadowing developed in |46|| and based on the light-cone 
Green function technique [47[. The propagation of a color dipole in a nuclear medium is described as motion in an 
absorptive potential, i.e. 



. dG{z2,r2;zi,ri) A^^G {z2,r2;zi,ri) 



dz2 



iya{l — a) 



krmnG{z2,r2;Zi,ri) 



iPA{z2,r2)(7gq{r2) 



G(z2,r2;zi,ri), (16) 



where the Green function G {z2,r2', zi,ri) describes the probability amplitude for the propagation of dipole state with 
size ri at the light-cone starting point zi to the dipole state with size r2 at the light-cone point Z2, and 

Q^a{l — a) +m'l 



2i/a(l - a) 
Then the shadowing correction to the amplitude has the form 
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Equation (J16p is quite complicated and in the general case may be solved only numerically |48j . However in some 
cases an analytic solution is possible. For example, in the limit of long coherence length, Ic 3> Ra, relevant for 
high-energy accelerators like the LHC, one can neglect the "kinetic" term ex Ar^G (z2, ''2; -^i, ^"1) in ([H]), and get the 
Green function in the "frozen" approximation |47| . 
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Then the shadowing correction (ITTI) simplifies to 
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If we neglect the real part of the amplitude and the longitudinal momentum transfer kmm (which is justified for 
asymptotically large s), and average over polarizations, then taking the integral over zi_2 "by parts" in (|19p . we get 
for the elastic amplitude 
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Another case where an analytical solution is possible is when the effective dipole sizes are small and the function 
<Jqq{r) may be approximated as 



,{r) « Cr^ 



(21) 



This approximation cannot be precise even at high virtualities Q^ and Mj, in DDVCS, since there are contributions 
of the aligned jet configurations mentioned above, which permit large dipoles even for large virtualities. Moreover, 
such aligned jet configurations of the dipole provide the main contribution to nuclear shadowing [431 ■ Nevertheless, 
for the sake of simplicity we use this approximation in order to estimate the magnitude of the shadowing corrections 
in the region xb & (10~^, 10""'^) . The approximation (|2ip is well justified on heavy nuclei. Namely, nuclear shadowing 
is independent of the form of the dipole cross section for large dipole sizes, above the saturation point, r^ > ^/Ql, 
where the typical value of the saturation momentum is Q^ ^ 1 GeV^ for heavy nuclei. Indeed, in this case the 



nucleus is "black". Therefore the shape of the dipole cross section matters only at r < 



0.16 /?7i. This size is 



sufficiently small for using the r approximation ((2T|) . Numerically, the approximation (|211) was tested in [48j-it was 



found that the discrepancy between the approximation (|21l) and the exact numerical solution of (|16p . changes the 
nuclear shadowing for DIS only within ten percent. We expect that within the same accuracy the approximation (|21|) 
is valid for DVCS. 

Then Eq. (J16p yields for W {z2, r2', zi,ri) the well-known evolution operator of harmonic oscillator, although with 
complex frequency 
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GLUON SHADOWING 



It has been known since [50, |51| that in addition to the quark shadowing inside nuclei there is also shadowing of 
gluons, which leads to attenuation of the gluon parton distributions. While nuclear shadowing of quarks is directly 
measured in DIS, the shadowing of gluons is poorly known from data |52l |53| . mainly due to the relatively large 
error bars in the nuclear structure functions and their weak dependence on the gluon distributions, which only comes 
via evolution. The theoretical predictions for gluon shadowing strongly depend on the implemented model-while for 
Xb ^ 10^2 they all predict that the gluon shadowing is small or absent, for xb ^ 10" 
wide range (see the review [53' and references therein). Some of the recent analysis [64' 

shadowing, that the unitarity bound |65| was severely broken. Since in this paper we also make predictions for the 
LHC energy range, the gluon shadowing corrections should be taken into account as well. 

The attenuation factor i?„ 



2 the predictions vary in a 
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where Gn {x,Q^,b) is the impact parameter dependent ghion PDF, was evaluated in the dipole approach in |54l.l55|. 
It was found convenient to evaluate Rg relating it to the shadowing corrections in DIS with longitudinally polarized 
photons, 
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where Act^ ^ = '^Z ~ ^'^1 ^ ^^ ^^^ shadowing correction at impact parameter b, and ct^ ^ {x, Q^) is the total 
photoabsorption cross section for a longitudinal photon. The process with longitudinal photons is chosen because the 
aligned jets configurations are suppressed by powers of Q^, so that the average size of the dipole is small, (r^) ~ 1/Q^, 
and nuclear shadowin g m ainly originates from gluons. 
As it was shown in [54l . |55| , 
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where pa{^^ z) is the nuclear density, and T (x^ Q^^ Az) is defined as 
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u — f cosh (riAz) + sinh (riAz) , 

w = (l + t^)sinh(f7Az) + 2tcosh(r2Az). 

Notice the importance of the Q-independent term in b^ , which controls the mean quark-gluon, or glue-glue separation 
at low scale. The magnitude of this term dictated by various experimental data ,63] and especially by data on 
diffraction 5J|, is rather large 6q ~ (0.65 6*6^)2. ^}y^^ leads to a small dipole sizes tq ^ 0.3fm and weak gluon 
shadowing. 

For heavy nuclei we may rely on the hard sphere approximation, pAii") ~ pyi(0)O {Ra — r)i and simplify (|24p to: 
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Acr-^'P (x, Q2) ^ j d%Aa'^'P{x,Q^,b) 



TTp^O) 

12 



n2R 

i dAzT{x,Q^,Az){lQR\-l2R\Az + Az 
Jq 



The results of the evaluation of the gluon shadowing are presented in Section IVIII 



VI. WAVE FUNCTIONS FROM THE INSTANTON VACUUM 

In this section we present briefly some details of the evaluation of the photon wavefunction in the instanton vacuum 
model (see [56l - l58| and references therein). The central object of the model is the effective action for the light quarks 
in the instanton vacuum, which in the leading order in N^ has the form [53, [53] 



S = 



d^x( — lnX + 2$2(2;) ~^(p + v~m~cLf(g><^- r,„ (g> fl) ?A j 



where Tm is one of the matrices, Tm = 1, ir, 75, ir75, ip and $ are the fields of constituent quarks and mesons 
respectively, N/V is the density of the instanton gas, v = w^7'^ is the external vector current corresponding to the 
photon, L is the gauge factor. 



L{x,z)^PexpUfdCv^{C)\ , 



(25) 



which provides the gauge covariance of the action, and f{p) is the Fourier transform of the zero-mode profile. 

In the leading order in N^.^ we have the same Feynman rules as in perturbative theory, but with a momentum- 
dependent quark mass fi{p) in the quark propagator 



S{p) = 



1 
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The mass of the constituent quark has a form 



Ai(p)=m + M/2(p), 



where ttt. « 5 MeV is the current quark mass, M « 350 MeV is the dynamical mass generated by the interaction with 
the instanton vacuum background. Due to the presence of the instantons the coupling of a vector current to a quark 
is also modified. 






ynom ^ _2Mf{p)^v^{q) + 0{q^) 
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Notice that for an arbitrary photon momentum q the expression for V^ depends on the choice of the path in ([^5]) 
and as a result one can find in the literature different expressions used for evaluations [59l462l|. In the limit p ^ 00 
the function f{p) falls off as ^ ^, so for large p ^ p~^, where p « (600 MeF)~^ is the mean instanton size, the mass 

of the quark /i(p) ~ m and the vector current interaction vertex V ~ v. However, we would like to emphasize that 
the wavefunction ^(/3, r) gets contributions from both the soft and the hard parts, so even in the large-Q limit the 
instanton vacuum function is different from the well-known perturbative result. 

We have to evaluate the wavefunctions associated with the following matrix elements: 
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where F is one of the matrices F = 7^, 7^75, cr^,y. In the leading order in Nc one can easily obtain 
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(29) 



The evaluation of ([29|) is quite tedious but straightforward. Details of this evaluation may be found in |59j . 

In what follows, we use Eqn (|29p for the evaluation of the initial and final wavefunctions ^m,/. Without any 
loss of generality, we may choose the frame with q = (A/^;, 0,0,0) for the time-like photon, and the frame where 
q = (0, 0, 0, Q) for the space-like photon, and after that make a boost to the Bjorken frame. However, we would like 
to emphasize that application of (P^ to the case of time-like photon should be done with care. Since the instanton 
vacuum model does not possess confinement, rigorously speaking Eq. (1291) may be applied only if the virtuality is 
below the quark pair production threshold, M^, ~ 4/i^(0) ^ 



(in case of pQCD wave functions the corresponding 
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Figure 1: DDVCS cross-section on the proton. Left: XB-dependence for different kinematical points. Right: i-dependence for 
different kinematical points. 



threshold is located at Mj, ^ Ami). Such behavior is characteristic to all models which have quarks as degrees of 
freedom but do not have built-in confinenient. 

The overlap of the initial and final photon wavefunctions in ([2|) was evaluated according to 



^f* (/3, r, -Af^) *(^) (/3, r, Q^) =^,1* (/?, r* , -Ml) h (/?, r, Q^) 
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where the summation is over the possible polarization states T — 7^,7^75,cr^i/- In the final state we should use 

, which is related to the reference frame with q'^ = 0, in which the components (j29p 



r„ = rn+n, 
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were evaluated. 
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VII. NUMERICAL RESULTS 



In this section we present the results of numerical calculations. While currently there is no data for DDVCS, we 
expect that similar to DVCS the experiments will be done in the region of large virtualities Q^. In that region we 
have Bjorken scaling, so all the model parameters such as basic cross-section uo and saturation radius Rq in Eq. ([9]). 
should depend on the Bjorken xb ■ A widely accepted parameterization which incorporates this feature is the GBW- 
type parameterization [27, 36-3^, |40] ■ The DDVCS cross-sections on the proton are shown in the Figure [T] Similar 
to DVCS, the DDVCS cross-section increases in the small- a: region like a;". As a function of t, the cross-section is 
exponentially decreasing. 

The Q^- and M^-dependence of the DDVCS cross-section on the proton is shown in the Figured As a function of 
the virtuality Q^, the DDVCS cross-section weakly depends on Q^ for Q^ < Mj,, but decreases like Q^^ for Q^ ^ M? 
Similar behavior is observed for M^da/dtdM^ for fixed Q^. Physically, such behavior is clear: the average dipole size 
(r^) is controlled by the wave function with largest virtuality, being (r^) ~ l/max((5^, M^). 

As one can see from the Figure [31 the shadowing correction is increasing towards small xb, and for xb ^ 10^^ 
the nuclear cross-section ratio decreases by a factor of two compared to the naive estimate da a ~ F'^(t)da]si- As a 
function of the momentum transfer t, the shadowing correction reveals a behavior qualitatively similar to the nuclear 
formfactor FA{t): it steeply drops at small-i and has zeros for some t. Notice, however, that the zero positions in 
the cross-section do not coincide with the zeros of the formfactor. This is a result of shadowing which suppresses the 
contribution of the central part of the nucleus and modifies the 6-dependence of the cross section compared to the 
formfactor. 
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Figure 2: DDVCS cross-section on the proton. Here we plot M^da/dtdM^ instead of da/dtdM^ in order to hide the trivial 
1/M^ in Eq. ([l]). Left: Q'^-dependence for different kinematical points. Right: A/^-dependence for different kinematical points. 
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Figure 3: The nucleus to nucleon cross section ratio for the DDVCS as function of different kinematical variables. Left: xb- 
dependence of the shadowing, t = tmin, Mj^ — 10 GeV . for different Q and A. From top to bottom: A = AQ and A — 208. 
Right: i-dependence, xb = 10"^, M^ = lOGeV^ 



The Q^- and M^-dependence of the shadowing is shown in the Figure|31 The Q^-dependence is similar to the DVCS 
case: the shadowing ratio is homogeneously increasing and for asymptotically large Q^ reaches 1. The shadowing is 
also increasing as a function of M^, but not so fast as a Q^-dependence. 

Concluding, we considered DDVCS on the proton and nuclear targets within the color dipole model. We found 
that the magnitude of the cross-section is small, of order a fe w p icobarns, and thus requires accelerators with high 
luminosities (e.g. future electron-ion colliders EIC and LHeC, [2l|,l22). The nuclear shadowing in the process is large 
and is important for analysis of DDVCS on nuclei. 
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